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flj ' The first physics results of the HARP experiment are presented. We emphasize the high 

I , performance of the forward part of the apparatus. The differential raw pion yield and its 

j~(. efficiency correction up to polar angles of 250 mrad are shown. The analysed setting is 

12.9 GeV/c incident protons in a 5% interacion legth aluminium target. 
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The HARP experiment'^ was designed to perform a systematic and precise study of hadron 
production for beam momenta between 1.5 and 15 GeV/c , for target nuclei ranging from 
hydrogen to lead. The detector was located at CERN, in the PS beam, and took 420 million 
events during the years 2001 and 2002. 

The physics goals of HARP are to make a measurement of the pion yield that will enable a 
quantitative design of the proton driver of a neutrino factory, and to improve the precision of 
atmospheric neutrino flux calculations. In addition, the energy-range is suitable to measure par- 
ticle yields for the prediction of neutrino fluxes for the MiniBooNE? - and K2K'^ experiments (8.9 
and 12.9 GeV/c respectively). To this end a collaboration was set up with these groups. Beam 
energy settings and dedicated targets were used to provide the most relevant measurements. 

A schematic layout of the apparatus is shown in Fig. ^ It is a large acceptance spec- 
trometer, with two distinct regions. A forward region (up to polar angles of about 250 mrad), 
where the main tracking devices are a set of drift chambers, and where particle identification is 
possible thanks to the combination of a threshold cerenkov, a time-of-flight wall and an electron 
identifier. In the large-angle region the main tracking and particle-id detector is a TPC, which 
is complemented by a set of RPC detectors for time of flight measurements. In addition, a com- 
plex beam instrumentation was set up in order to establish the nature of the incoming particle, 
given the rather impure beam. The main subsystems are three small timing detectors and two 



cerenkovs. The beam time detectors are also used to estimate the interaction time at the target, 
located inside the TPC. 
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Figure 1: Schematic layout of the HARP spectrometer. 

Given the immediate interest of the MiniBooNE and K2K experiments in our results, the 
HARP Collaboration has decided to start the data analysis with these two attractive items. In 
particular, we present in this article a first analysis of the data dedicated to the K2K measure- 
ment, using for this purpose the forward region of the detector, which is able to cover by itself 
the K2K requirements. 

The organisation is as follows. Section[21introduces the motivation of this analysis. In section 
I21the tracking capabilities of the forward spectrometer are presented. Section|3]is devoted to the 
particle identification subsystems and their performance. Finally, the data analysis is presented 
in section 13 and the conclusions in sectional 

2 Motivation of this analysis 

One of the main systematic errors on the neutrino oscillation parameters measured by the K2K 
experiment comes from the uncertainty on the far/near neutrino flux ratio. This ratio depends 
on the differential pion production cross section, which is essentially measured by a threshold 
cerenkov (pion monitor) above E^ = 1 GeV, but must rely on Monte Carlo simulations below 
that energy. Unfortunately, the oscillation peak is located at Ei, ~ 0.55 GeV. HARP could 
precisely estimate the pion differential cross section at these energies and subsequently reduce 
that systematic error. 

Fig. 13-left shows the (p, 0) distribution for pions producing neutrinos in the bin containing 
the oscillation peak (0.5 — 0.75 GeV). The relevant phase space is l<p<8 GeV and 9<250 mrad, 
which is fully covered by HARP's forward region as shown in Fig. |2J 

3.5 million useful events were taken by HARP with incident 12.9 GeV/c protons and an 
exact replica of the K2K target, which is a 80 x 3 cm aluminum tube (2 interaction lengths). 
Another 6M events were collected with a 5%A aluminium target ("K2K thin target") and a 
similar beam in order to decouple reinteraction and absorption effects from pure production. 
The data analysis presented in this article is based on this particular setting. 



3 Forward Tracking 

Tracking of forward going particles is done by a set of drift chambers (NDC) placed upstream 
and downstream of the dipole magnet (DIP). The NDC-DIP combination allows the momentum 
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Figure 2: On the left, {p,0) distribution for pions producing neutrinos in the bin containing the oscillation peak 
(0.5 — 0.75 GeV). This figure is a courtesy of the K2K Collaboration. It has been obtained with a full simulation 
of the K2K pion production and decay chains. Center and right, geometrical acceptance of HARP's dipole 
magnet for particles arriving to the first NDC module (NDCl). On the center, P dependence for particles with 
\6x\ < 200 mrad. On the right, 9^ dependence for particles with P > 1 GeV/c. 



measurement by matching track segments located at either side of the magnet. 

The chambers were recuperated from the NOMAD experiment and their properties have 
been described elsewhere'^'. Each NDC module contains 4 chambers, and each chamber 3 planes 
of wires with tilted angles —5°, 0° and 5°. The single wire efficiency is of the order of 80%, and 
the spatial resolution approximately 340 fim. " 

The reconstruction algorithm builds 2D and 3-D track segments in each NDC module (12 
hits maximum), which are fitted to a straight line model via a Kalman Filter fit'^. Afterwards, 
all the possible matching combinations of tracking objects (including unused hits) belonging to 
different modules are performed in order to obtain longer tracks. 

The momentum measurement is done associating 3-D downstream segments with 3-D and 
2D segments in the upstream module (NDCl). This asymmetry is needed to compensate the 
low tracking efficiency of the upstream chambers, due to the presence of a unique module and to 
the higher hit density (proximity to the target), which provokes pattern recognition confusion. 
Fig. ISJ-bottom shows the momentum and angular resolutions for 3D-3D up-down matches. 

The small track separation in NDCl induces a huge correlation between particles that implies 
a hadron model dependent tracking efficiency. This is a potential source of systematic error as 
shown in Fig. |31-top. The average efficiency is of the order of 65% when one consideers only 
3.D-3-D up-down matches. The situation improves when 2D-3D matches are included. Recent 
studies with a third type of tracks, built by matching a 3-D downstream track with the vertex, 
show a considerable efficiency recovery. As these tracks are independent of NDCl, the total 
efhciency (3-D-3-D + 2D-3D + vertex-3D) is nearly model independent. 



4 Particle Identification 



Particle identification in HARP's forward region uses the information from the time of flight 
system (TOF), the threshold cerenkov (CKOV) and the electron identifier (EID). Pion/proton 
separation is provided by TOF up to 4.5 GeV/c, and by the CKOV above 3 GeV/c. Elec- 



"The chambers had a superior performance in NOMAD (95% hit efficiency and a resolution a factor of two 
better). This is mainly due to the use of different non-flammable gas mixture (y4r(90%) — C02(9%) — (77^4(1%)), 
and volatage settings ( sense wires held at + 1300V and potential wires at -2900V). However, one should stress 
that this spatial resolution is more than sufficient for HARP physics 



tron/pion separation is covered by the CKOV below 3 GeV/c and by the EID above 2 GeV/c. 
Finahy the kaon contamination can be estimated with the CKOV above 3 GeV/c and with the 
(TOF) below this energy. 
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Figure 3: On the top, tracking efficiency as a function of momentum (left) and angle (right). The bands rep- 
resent the maximum difference between 3 Monte Carlo hadron generators (MCI, MC2, MC3). The lower band 
corresponds to 3D-3D up-down matches while the upper band includes also 2D-3D up-down matches. On the 
bottom- left, momentum resolution for data (dark) and Monte Carlo (light). On the bottom-right, angular reso- 
lution for Monte Carlo. Both correspond to 3D-3D matches, with no vertex constraint included. 



4-1 Time of Flight System 

Particle identification by time of flight in HARP's forward region relies on the combination of 
particle momenta {p) and track length (L), measured with the forward spectrometer, and the 
time-of-flight between a start signal (to) from the beam time detectors and a stop signal {tyj) 
from the TOF wall, placed at about 10 meters downstream from the target. The mass of a 
particle can be computed from these quantities, m? = p^ ■ [((t^ — to) • c/L)^ — 1]. 

The current TOF wall and to resolutions are 150 and 70 ps respectively. However, the to 
resolution used for the analysis shown in this paper was above 200 ps, leading to a combined 
TOF resolution of about 270 ps, which is already better that the design value of 300 ps. 

The TOF is able to provide a vr/p separation of ~ Sc at 3 GeV jc, as can be seen in Fig. 
I^left. In practice, the TOF system will provide reasonable vr/p separation up to 4.5 GeV jc. The 
tt/A; capabilities of the detector at low energies (<3 GeV jc) are being studied at the moment^. 



''It has been recently proposed to use unbiased golden events (all beam time detectors firing to gain in io 
resolution) and tracks (crossing the overlap region between counters in the TOF wall to increase the tu, resolution) 
in order to estimate the vr/fc ratio below 3 GeV jc. The Iq resolution obtained in this case is close to 130 ps. Anyway, 
the kaon contamination at these energies is known to be very small. 
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Figure 4: On the left, mass squared as measured by the TOF system for 3 GeV jc beam particles. On the center, 
photoelectron yield in the cerenkov produced by beam pions, electrons and protons in a 5 GeV jc no target run. 
On the right, hadron and electron yields as a function of E\/E and E/p for beam particles in a 3 GeV jc no target 
run. E is the sum of the energies deposited in the two calorimeter modules (E\ + E2) and p is the momentum. 
In all cases pure samples are selected using the beam detectors. 



^.2 Threshold cerenkov 

The cerenkov detector is able to identify electrons below the vr saturation (~ 4 GeV/c). Pions 
and electrons can be separated from protons and kaons below the kaon threshold at ~ 9 GeV/c, 
and protons can be discriminated above this energy. Fig. [IJcenter shows the photoelectron yield 
for beam pions, electrons and protons in a 5 GeV/c no target run. At 5 GeV/c pions are close 
to the saturation regime and cannot be distinguished from electrons, which on the other hand 
are very rare at these energies. However, ir/p separation is highly pure and efficient at 5 GeV/c 
(and above). 

4.3 Electron identifier 

The electron identifier is made of two calorimeter planes reused from the CHORUS experiment 
and described elsewhere'^. It was designed to provide electron-pion separation when low energy 
(< 3 GeV/c) charged pions, accompanied by knock-on electrons, are occasionally identified as 
electrons by the cerenkov counter. It also serves to identify electrons at high energy, when 
the cerenkov has lost its pion/electron separation capabilities. The detector performance is 
summarised in Fig. 01-right. 

5 The analysis 

One sixth ( 1 million events ) of the "K2K thin target" data has been analysed. The unnormalised 
pion production differential cross section can be computed as follows'^: 
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where the index i corresponds to {p, 9) bins, e'^'^^ is the geometrical acceptance, eY"''^^ is the 
tracking efficiency, ej is the pion identification efficiency, jyf is the pion purity and N^ is the 
observed pion yield. The pion purity is defined as r]J = {N[ — N^ ^)/N'^ 



where A'^^ ^ is pion 

misidentification background. Fig. El shows the raw pion yield corrected by the tracking and 
pion identification efficiencies ( tratk tv ^i) as a function of p and 9. 



"Bin migration effects axe not considered for the moment. Preliminary studies show a small effect. 



The current activities in the context of this forward analysis follow the line of recovering 
tracking efficiency by using the matching between the downstream tracks and the vertex. The 
most recent results obtained with the improved tracking are very encouraging as they are nearly 
model independent. This will allow a considerable reduction of the systematic error. 
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Figure 5; Raw pion yield corrected by tracking and pion identification efficiencies. On the left, the p/n misiden- 
tification background for the TOF system is also shown (for the three hadron generators). The number of Monte 
Carlo events used to compute the tracking efficiency was of the order of 50000. This explains the large statistical 

errors in these figures. 



6 Conclusions 

We have described the current performance of the HARP apparatus and our first physics analysis. 
This analysis is based on specific data taken by HARP to improve the neutrino flux calculation 
for the K2K experiment. Particle identiflcation and tracking of forward going particles is well 
under control. This has permitted an initial estimation of the differential raw pion yield and its 
efficiency correction. The analysis is still in progress and will lead to conclusive results shortly. 
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